Dynamical mechanisms of DC current generation in driven Hamiltonian systems 
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Recent symmetry considerations (Phys. Rev. Lett. 84 2358 (2000)) have shown that dc currents 
may be generated in the stochastic layer of a system describing the motion of a particle in a one- 
dimensional potential in the presence of an ac time-periodic drive. In this paper we explain the 
dynamical origin of this current. We show that the dc current is induced by the presence and 
desymmetrization of ballistic channels inside the stochastic layer. The existence of these channels 
is due to resonance islands with non-zero winding numbers. The characterization of the flights 
dynamics inside ballistic channels is described by distribution functions. We obtain these distribution 
functions numerically and find very good agreement with simulation data. 
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Transport in driven systems has received widespread 
interest for several years because of its potential appli- 
cability to nonequilibrium processes fl. One canoni- 
cal model reduces to the motion of a particle in a one- 
dimensional space-periodic potential in the presence of 
friction and a time-dependent stochastic force x(t) ■ 11 
x(t) contains correlations, a nonzero current may be re- 
alized even in the case of zero average < x(t) >— |Q. 
Despite an enormous accumulation of results in this area 
H we are still lacking a full understanding of the mi- 
croscopic mechanisms of current occurrence. If such an 
understanding is realizable, it should make use of the true 
dynamical evolution of the system rather than of prop- 
erties of equations for probabilities. A first step in such 
a direction requires to separate the essential time corre- 
lations from the pure Gaussian white noise in The 
simplest way is to assume that x(t) = E(t) + £(t) where 
E(t) = E(t + T) is a time-periodic function with zero 
mean < E(t) >— and is a Gaussian white noise 
term. 

The next step is to skip the £(t) term which leaves us 
with a regular dynamical problem. In || such a case was 
considered and the relevant space-time symmetries of the 
dynamical problem have been obtained. It was shown 
that a breaking of those symmetries leads to a nonzero 
dc current. The mechanism of current occurrence for 
the dissipative case was identified with a desymmetriza- 
tion of attractor basins. In contrast the nondissipative 
(Hamiltonian) case is much less understood Q]. A time- 
dependent Hamiltonian system is usually nonintegrable 
pj . A strong dc current component was found in the cor- 
responding stochastic layer of such a system pj . While its 
presence or absence was clearly connected to the above 
mentioned absence or presence of symmetries, the dy- 
namical nature of directed transport in the stochastic 
layer lacks full understanding. The importance of this un- 
derstanding can be seen from e.g. results in || where ki- 
netic equations for probability functions have been stud- 
ied. In particular it was found that the approaching of 



the Hamiltonian (dissipationless) limit leads to an in- 
crease of the dc current value by 2-3 orders of magni- 
tude. Thus the description of the dynamical mechanisms 
of directed current generation in the stochastic layer of a 
driven Hamiltonian system will provide with very useful 
information for dissipative systems as well. 

Let us consider the canonical example of a parti- 
cle moving in a spatially periodic nonlinear potential 
U(x) = — cosx under the influence of time-periodic zero- 
mean force E(t). The Hamiltonian and the equation of 
motion are given by: 



H = 



P 



cos a; 



xE(t) 



SIM 



+ E(t) . (1) 



Here p and x are canonically conjugated momentum and 
coordinate and x = d 2 x/dt 2 . 

We restrict our consideration to the choice 



E(t) = E x cos(i) + E 2 cos(2i - 



(2) 



According to || for E<x ^ and ^ 0, n all possible 
symmetries which yield zero dc current are broken. Note 
that the phase space dimension d of (Q) is d = 3. 

In the case of a nonzero field E(t) the phase space of 
(|l|) is characterized by the presence of a stochastic layer 
which originates from the destroyed separatrix of the un- 
driven system For <j> — 0,7r this layer is invariant 
under the transformation (ja — » — p, t — > — t, x — > x). 
At the same time the average velocity for any trajectory 
in this layer vanishes, so we find zero dc current. The 
symmetry will be broken when tuning cj) away from the 
values (0,7r). The stochastic layer will deform. Most im- 
portantly any trajectory in the layer will then be charac- 
terized by a nonzero value of the average velocity. Due to 
ergodicity inside the layer this value will be unique for all 
trajectories from the layer. While the fact that it may 
become nonzero is understandable using the symmetry 
analysis, its appearance and magnitude is due to dynam- 
ical mechanisms of motion inside the stochastic layer. In 
this paper, we show that the dc current is induced by 
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the presence and desymmetrization of ballistic channels 
inside the stochastic layer. The existence of these chan- 
nels is due to resonances. The characterization of the 
realization of flights inside ballistic channels is described 
by distribution functions. We obtain these distribution 
functions numerically and find very good agreement with 
simulation data. 

System (Q) has a mixed phase space, which contains 
of chaotic areas and regular resonance islands . These 
islands are impermeable for chaotic trajectories and, at 
a first glance, may be excluded from the consideration 
of the phase space flow inside the stochastic layer. In 
reality the phase space topology inside the stochastic 
layer is very complex precisely at the boundary between 
chaotic and regular regions M. Close to resonances the 
stochastic layer shows up with a hierarchical set of can- 
tori, which form partial barriers for a trajectory from the 
layer. Due to the presence of these barriers a chaotic 
trajectory can be trapped for a long time near a reg- 
ular island resonance. This trapping or sticking effect 
leads to the appearance of strongly nonergodic episodes 
during the overall chaotic motion. Regular islands are 
characterized by a corresponding rational winding num- 
ber to = Ax/T which defines the distance Ax travelled 
during one period T = 2n of the drive. If the winding 
number lu is nonzero, the corresponding sticking episode 
of the chaotic trajectory is a ballistic-like unidirectional 
flight. For u = the sticking episode corresponds to 
trapped oscillations. 

Thus the complicated evolution of a trajectory in the 
stochastic layer can be subdivided into several parts [Q . 
The first one is a fast diffusion in the bulk of the layer, 
while the other ones are stickings to the above mentioned 
regular islands and correspond to propagation in ballistic 
channels. The switching from the diffusion process into a 
ballistic flight will be described by some probability dis- 
tribution. The same will be true for the actual residence 
or sticking time inside a given channel. We will show for 
the cases studied that the fast diffusion alone is not capa- 
ble of explaining the observed dc current. The main point 
is, that the leading mechanism of current generation in 
the stochastic layer is related to the desymmetrization of 
the strongly non- stochastic part of the overall stochas- 
tic dynamics inside the layer. Note that our kinetic en- 
ergy choice p 2 /2 in (|l]) implies that the stochastic layer is 
bounded in p, so we will always expect ballistic channels 
to appear. 

Let us study the case of weak driving E\ — 0.252 and 
Ei = 0.052. A Poincare map of the phase space flow for 
<p = is shown in Fig. 1(a) M. The main stochastic layer 
(central location) shows up with zero average velocity due 
to symmetry arguments (Fig. 2). The large hole in the 
middle of this layer corresponds to regular trapped mo- 
tion in the well of U(x). Additional resonances are seen 
above and below the central layer. These thin ballistic- 
like but yet stochastic channels have no overlap with the 



central layer. 
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FIG. 1. Poincare map for (a) <f> — and (b) <j> = 7r/2. 

A weak asymmetry <p = tt/5 leads to a slight deforma- 
tion of the main stochastic layer and to a desymmetriza- 
tion of the overlap of the chaotic layer with higher-order 
resonances and to the appearance of a positive current 
in the system. Note that it still does not overlap with 
the thin ballistic channels seen in Fig. 1(a). Most impor- 
tantly we observe a nonzero average velocity < x >~ 0.05 
(Fig.2). 

Further increase of the asymmetry, <p = 7r/2, results 
in an overlapping of the main stochastic layer with the 
upper ballistic resonance Fig. 1(b). Note that at the same 
time the lower ballistic resonance is not overlapping. The 
average velocity increases to < x >« 0.2, which is four 
times larger than the result for cf) = ir/5. Standard har- 
monic mixing theories (see 0, ||) would predict a de- 
pendence < x >^ sin0 and thus only an increase by a 
factor of 1.7. 

In the x(t) curves in Fig.2 we observe many ballistic 
flights. For one of them an inset shows the corresponding 
Poincare map result, which verifies that these flights cor- 
respond to stickings of the chaotic trajectory to the upper 
balistic resonances. A zooming of the x(t) curves shows 
self-similarity, i.e. the seemingly random dynamics be- 
tween observable long flights is actually again composed 
of shorter flights and seemingly random dynamics etc (cf. 
insets in Fig.2). 

In order to quantify our analysis of the symmetry bro- 
ken dynamics we compute the distribution of travelling 
times of 'uniform' flights to the left P_(tf) and to the 
right P + (tf) separately. Here 'uniform' means no change 
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of direction of motion [|10| . For each separate flight we 
note both the time tt spent in this motion and the dis- 
tance x f travelled. The dependence of Xf on tf is shown 
in the inset of Fig. 3(a) for (j> = ir/2. Similar to the other 
cases c/> = 0,7r/5 we observe a simple fork-like structure. 
This is due to the fact that any considerable distance 
covering in the stochastic layer is realized through flights 
while sticking to the boundary of the stochastic layer. 
The slopes in the inset of Fig. 3 (a) are given by the corre- 
sponding winding numbers of the layer boundaries. Note 
that the two fork parts merge at values of tf « 10T. In 
principle other ballistic channels with different winding 
numbers might be present. Here they were too weak to 
be detected. 
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FIG. 2. x(t) for cj> — 0,7r/5,7r/2 (lower, middle and upper 
curves respectively). Left upper inset: Poincare map for bal- 
listic flight with <f> = 7r/2 as indicated by arrow. Right inset: 
zoom of x(t) for (j> = t/2. 

In Fig. 3 (a) we show the corresponding distribution 
functions P±(tf) (again for <f> = tt/2). They are obtained 
by counting the number of flights with tt falling into a 
time window of size T. For tt < 10T we observe expo- 
nential dependence of P± on tt . These short flight distri- 
butions are in fact independent of the direction of flight. 
Skipping all longer flights would lead to the prediction 
of nearly zero average velocity (restricting to flights of 
length tf < 10T yields about one percent of the numeri- 
cally observed current). Thus the desymmctrization will 
manifest itself for longer flights. For tt > 10T a crossover 
to a power law P± ~ t?* takes place. Here we find a 
significant desymmetrization for <fi — 7r/5,7r/2. Estimat- 
ing the exponents [jll| we find for <j> = 7r/5: a_ « 2.5, 
a + 2.4 and for <j> = tt/2: a_ « 3.7, a + ps 2.3. It 
is worthwhile noting that a < 3 implies unidirectional 
anomalous diffusion with diverging second moments of 
P(t). The flights are coined Levy flights in such a case 



Following the continuous-time random walk (CTRW) 
formalism Jl3| we propose a generalized asymmetrical 
flight model capable of reproducing the above results. 
The applicability of the CTRW model follows from the 
assumption that the presence of a random phase with 
fast decaying correlations leads to the absence of corre- 
lations between consecutive flights, since they are almost 
always separated by dispersive chaotic motion. 
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FIG. 3. P+(t f ) (solid line) and P-(t f ) (thick dashed line). 
Insets: Xf versus tf. For parameters see text. 

Assume that there exist N different resonances with 
winding numbers Wi, i = 1,...,N. Every resonance 
is characterized by a probability distribution function 
(PDF) of sticking time Si(t). After finishing a random 
phase event, the probability of sticking to the ith reso- 
nance is pi, Yli=i Pi = 1- The random phase residing 
time is characterized by a PDF S r (t). All functions Si(t) 
and S r (t) must have finite first moments, due to the Kac 
theorem about finitcness of recurrence times in Hamilto- 
nian systems 0. With these definitions we obtain the 
following expression for the current: 



TtLiUiPijtj) 

<*r> 



J 



(3) 



where (t t ) = f tPi(t)At. 

For the above discussed cases of <j> = 7r/5,7r/2 we find 
only two relevant ballistic channels - one with positive 
winding number and a second one with negative wind- 
ing number. In order to properly obtain S(t), we note 
that our numerically obtained function P(t) consists of a 
lot of short 'flights' as defined through the numerics p0[ . 
These may be cither stickings to islands with zero wind- 
ing number or chaotic motion. We observe that for flight 
times tf > 10T only ballistic flights with nonzero wind- 
ing number are obtained. So the functions S± (t) may be 
easily obtained from P±(tf) by cutting the central part 
t < 10T out and properly normalizing. In this case the 
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expression for the average current simplifies to 



J = 



1 



«(1 + /) 



(w + (t+) + /w_(t_)) 



(4) 



where the two constants k and / can be obtained from 
the total time of a simulation T to t and the numbers N± 
of ballistic flights, k = T tot /(N++N-) and / = A_/A+. 

For <p — 7r/5 we obtain from the numerical runs / ps 
0.57, k ps 1900, (*+) ps (i_) ps 220 and cj+ = 10/6 ps 1.67, 
u}— = — 1.5. In this case of weafc desymmetrization the 
main source of a nonzero current is the different proba- 
bility to enter a right or left going flight because / ^ 1. 
At the same time the average flight times in both ballis- 
tic channels are nearly identical. With the help of (|J) we 
find J ps 0.056 which is close to the numerically observed 
value 0.05. 

For the case 4> = vr/2 wc find / ps 0.16, K ps 2600, 
(t+) ps 400, ps 150 and uj + = 2, w_ ps -1.4. Note 
that the above discussed overlap with the upper reso- 
nance yields a further strong desymmetrization in the 
probability to realize a left or right going flight, and in 
addition the average flight times in both channels signif- 
icantly differ. Expression (Q) yields J ps 0.22 which is in 
good agreement with the numerically observed value 0.2. 

For stronger driving amplitude E\ = 3.26, E2 = 1.2 
and <fi — 7r/2 we obtain an average velocity (x) ps 0.85. 
The corresponding Xf(tf) dependence and the PDFs 
P±(tf) are shown in Fig. 3(b). The Xf(tf) dependence 
shows that more than two ballistic channels are involved. 
The asymmetry of the PDFs at short flight times in- 
dicates that a considerable number of left-going flights 
becomes dominating at short times, in agreement with 
the tendency of the previous results. This makes the 
application of the simplified sum rule (^) impossible, in- 
stead the original definition (||) should be used. Care- 
ful analysis of the structure of the stochastic layer shows 
that relevant resonances become embedded in the bulk of 
the stochastic layer. While these structures are of rather 
small size, they are frequently visited. A restriction to 
short flights tf < 10T now yields a considerable nonzero 
current which is however negative, i.e. opposite to the 
total current value. Again the long ballistic flights are 
necessary in order to properly obtain the observed cur- 
rent value. 

In summary, we have explained the dynamical mecha- 
nisms of current appearance in driven Hamiltonian sys- 
tems inside the stochastic layer with broken time-reversal 
symmetry. The key source of such a directed transport 
is the desymmetrization of flight probabilities in ballis- 
tic channels inside the layer. As it follows from our sum 
rule ^ the resulting current depends among other pa- 
rameters on the average time spent in a ballistic channel. 
These times will sensitively depend on control parame- 
ters of the system if the exponent a becomes less than 3. 
In such a case small changes may significantly alter the 
current value as shown above. 



A recently proposed geometric approach of counting 
areas and winding numbers is in principle also capable 
of obtaining the observed mean value for the current 
]l4| . This approach may also require sophisticated stud- 
ies of the fractal structure of the chaotic layer. It rep- 
resents a nontrivial complementary result, since it, al- 
though not explaining the dynamical mechanisms of cur- 
rent generation, is capable of obtaining the average cur- 
rent value (provided the sums in Eq.(3) of (l4| converge 
fast enough). 

We thank M. Fistul and A. A. Ovchinnikov for useful 
discussions. 
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